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ABSTRACT 

We extend a chemical evolution model relating galaxy stellar mass and gas-phase oxy- 
gen abundance (the mass-metallicity relation) to explicitly consider the mass-dependence of 
galaxy gas fractions and outflows. Using empirically derived scalings of galaxy mass with 
halo virial velocity in conjunction with the most recent observations of z ^ total galaxy 
cold gas fractions and the mass-metallicity relation, we place stringent global constraints on 
the magnitude and scaling of the efficiency with which star forming galaxies expel metals. We 
demonstrate that under the assumptions that metal accretion is negligible and the stellar initial 
mass function does not vary, efficient outflows are required to reproduce the mass-metallicity 
relation; without winds, gas-to-stellar mass ratios > 0.3 dex higher than observed are needed. 
Moreover, z = gas fractions are low enough that while they have some effect on the magni- 
tude of outflows required, the slope of the gas fraction-stellar mass relation does not strongly 
affect our conclusions on how the wind efficiencies must scale with galaxy mass. Because 
theoretical descriptions of the mass loading factor ?]„ = M^^/AIsfr, where M^ is the mass 
outflow rate and A/sfr is the star formation rate, are often cast in terms of the depth of the 
galaxy potential well, which is in turn linked to the host halo virial velocity Uvir, we use one 
of the latest abundance matching analyses to describe outflow efficiencies in terms of v^ir 
rather than stellar mass. Despite systematic uncertainties in the normalization and slope of 
the mass-metallicity relation, we show that the metal expulsion efficiency Cw = (^w/^g)?7w 
(where Z^ is the wind metallicitiy and Zg is the interstellar medium metallicity) must be both 
high and scale steeply with mass. Specifically, we show that Cw ^ 1 and (^ oc v~^'^ or steeper. 
In contrast, momentum- or energy-driven outflow models suggest that ri^j should scale as u^J 
or v^^^, respectively, implying that the Z^-AI^, relation should be shallower than the Zg-M^ 
relation. 

Key words: ISM: abundances — ISM: jets and outflows — galaxies: abundances — galaxies: 
evolution — galaxies: fundamental parameters — galaxies: ISM 



1 INTRODUCTION 

Star- forming galaxies follow a tight (~ 0.1 dex scatter) correlation 
between their gas phase oxygen abundance (hereafter referred to as 
"metallicity") and stellar mass JTremonti et al.ll2004r) . This mass- 
metallicity relation is primarily understood to be a sequence of 
oxyg en suppression, rather than enrichment rTremo nti et alj|2004l ; 
lOalc anton 2007; Erb 2008; Finlator & Dave 2008). The produc- 
tion of oxygen traces the production of stars, implying that the ob- 
served trend in the oxygen-to-gas ratio reflects either a trend in the 
galaxy gas-to-stellar mass ratio or in processes that affect gas-phase 
metals but not stars. A consensus is emerging that although galaxy 
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gas fractions can and do affect the mass-metallicity relation, if the 
stellar initial mass function (IMF) is the same in all galaxies, then 
outflows that are more efficient at removing metals from low-mass 
galaxies are required in order to reproduce the observation s (e.g., 
ll3alcantoiJl2007l : lFinlator & Davdl2008l : ISDitoni et al.ll2010l) . How- 
ever, the global properties of these outflows and the physics under- 
lying how star formation drives them are not well understood — and 
winds are expensive and difficult to observe directly. In this paper, 
we incorporate the most recent observations of galaxy gas fractions 
and the mass-metallicity relation at z ~ into a simple chemical 
evolution model to explore what constraints can be placed on how 
the efficiencies and composition of star formation driven galactic 
winds scale with galaxy stellar mass and halo virial velocity. 

Several analytic studies have concluded that star formation 
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driven outflows are cru c ial to reproducing the observed mass- 
metallicity relation. lErbl ( 120081) used a simple analytic chemical 
evolution model to argue that the star formation rate, AfsFR, and 
the outflow rate, Af„, should be roughly equal. While Mw and 
the gas accretion rate Mace vary with the star formation rate (and 
thus gas fraction), rj^ = M„/Msfr and 77a = Macc/MsFR are 
constant universal parameters, a common p ractice in analytic m od- 
els of galaxy chemical evolution (see also ISamui et al.ll2008l , and 
references therein). Though models specifically aimed at duplicat- 
ing observa tions of the mass-m etallicity relation commonly assume 
Z^ = Zg. lPalcantonl ( 120071) argues that metal-enriched outflows 
(those comprised predominantly of Type II supernova ejecta, and 
thus with Zw > Z^) are r equired if the rate of gas accretion is to be 
reasonable. More recentlv. lSpitoni et alj ( l2010l) have argued that the 
z = mass-metallicity relation together with gas fractions derived 
by inv erting the Kennicutt-Schmidt (K-S. lKennicuttll998l : ISchmidll 
Il959h law imply that not only are outflows required, but that they 
must be more efficient at rem oving metals from low-m ass galaxies 
than from more massive ones. lFinlator & Davg 1 120081) drew a sim- 
ilar conclusion by analyzing a suite of cosmological smoothed par- 
ticle hydrodynamics (SPH) simulations evolved with Gadget-2 
dSpri ngel 2005) in conjunction with detailed analytic models. They 
showed that, in general, Zg oc r?^" ^ for t^w 3> 1. Their favored 
model that reproduces the lErb et al.l 2 ~ 2.2 mass-metallicity re- 
lation is one in which r]^ oc a~^, where a is the galaxy velocity 
dispersion|j In this simulation, a~^ oc Afj^j^ oc M^ , which 
naturally explains why this t^w scaling is able to reproduce a mass- 
metallicity relation with Zg oc M^''*. These simple scaling rela- 
tions highlight a link between a galaxy's stellar mass, its halo mass, 
and its potential well: wind models aimed at successfully reproduc- 
ing the mass-metallicity relation also need to correctly reproduce 
(or incorporate) the A/i-A/haio relation. Moreover, this analysis 
shows that the mass-metallicity relation is a potentially powerful 
tool for constraining how star formation driven outflows scale with 
galaxy and halo properties; this is particularly interesting as such 
scalings are currently not well constrained through either direct ob- 
servations or theoretical considerstions. 

On the other hand, several models focus instead on the effi- 
ciency of star formation as a function of stellar mass. In such mod- 
els, an increase in the star formation efficiency with galaxy mass — 
without the need for outflows — is sufficient to reproduce the ob- 
served mass-metallicity relation JCalura et alJl2009f ). lBrooks et al.l 
1 20071) used a set of SPH simulations evolved with Gasoline 



I Wadslev et al.l 12004 ) — a nd therefore a di f ferent recipe for star- 
formation feedbaclQ than iFinlator & Davg ( 120081) — to argue that 
preferentially expelling gas from the low-mass galaxies is insuffi- 
cient for reproducing the observed mass-metallicity relation. These 
authors claim that it is instead the reduced star-formation efficiency 
(and thus differences in galaxy gas fractions) induced by such feed- 
back that is primarily responsible for driving the relation's mor- 



^ This parameterization is motivated by the observations of | MaitiiJ 120051) 
and the theory of momentum-driven winds jMuiTav et alJIlOOa) ; see ^2.3 
for more details. 

^ Because of the resolution of cosmological SPH simulations, star- 
formation feedback must be included using "r ecipes" instead o f directly 
modelling the underlying physics. The winds in lFinlator & Dava s simula- 
tions are implement ed by physica lly moving gas particles away from star- 
forming regions. In iBrooks et al.r s simulations, star formation thermally 
heats neighboring pailicles. In both prescriptions, the relevant particles ai'e 
not a llowed to inte ract hydrodynamically (Finlator & Dave) or radiatively 
cool IBrooks et alj) for some physically-motivated amount of time. 



phology. In the context of the mass-metallicity relation, variations 
in the star formation efficiency affect galaxy gas fractions (as well 
as the Af*-A'fhaio relation). We do not directly address star for- 
mation efficiency here because we take both galaxy gas fractions 
and the A'/*-Afhaio relation as givens rather than something to be 
constrained by the model; we discuss in Appendix iBl the implica- 
tions our choices for these relations have on how star formation 
efficiency varies with galaxy mass. 

Finally, letting the IMF (and thus the amount of oxygen 
produced per unit stellar mass) vary with galaxy mass pro- 
vides a straightforward way to reproduce t he mass-metallicit; 



relation JTinsley 1974; Koppen et al. 20071; iRecchi et 
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ICalura & Mencill2009l ; ISpitoni et alj|2010h . We primarily assume 
here that the IMF is the same in all galaxies, with a brief discussion 
in ij |4.2| of how uncertainties in yields in the presence of a variable 
IMF affect our results. In general, if the IMF is top-light in low- 
mass galaxies then this will imply that outflow efficiencies do not 
need to scale as steeply with mass as suggested by the non- varying 
IMF case. 

We apply here a simple model with which to understand the 
mass-metallicity relation to the mass-metallicity relation at 2 ~ 0, 
where external constraints such as gas fractions and the stellar mass 
function are best measured. We do not assume a particular form 
for the mass-metallicity relation; we instead base our conclusions 
on the range of parameter space allowed by the range of system- 
atic uncertainties in interpretting strong nebular emission lines as 
oxygen abundances (ij |2.1b . Our main simplifying assumptions are 
that the metallicity of gas accreted from the intergalactic medium is 
negligible and that the nucleosynthetic yield is constant with galaxy 
mass. With these constraints and assumptions, the only free param- 
eters are those describing outflows, which we are able to describe 
as a function of halo virial velocity. 

This paper is organized as follows. In §[2] we discuss the 
relevant observations. The slope and normalization of the mass- 
metallicity relation strongly affect the interpretted properties of 
galaxy winds. Unfortunately, while there exist exquisite data on 
emission line ratios of star forming galaxies at 2; = 0, the correct 
way to interpret these line ratios in terms of oxygen abundances 
is not agreed upon; we therefore consider several measurements 
of the z = Q mass-metallicity relation, as outlined in i) |2.1| It is 
commonly assumed in chemical evolution models, and we assume 
here, that the gas is well-mixed; we address the differences be- 
tween galaxies' cold gas resevoirs and the gas traced by star for- 
mation in §[2J2](see also AppendixlAJl. As the purpose of this paper 
is to place constraints on how galaxy outflows scale with galaxy 
mass, we briefly outline observed properties of galaxy outflows 
(and theoretical models thereof) in §|23] We lay out the formal- 
ism in i]3.1| and its derivation in Appendix |C] along with how we 
connect galaxy stellar masses to host halo properties (ij |3.2b . In §|4] 
we show how gas dilution and outflows must combine in order to 
yield the observed mass-metallicity relation, and what this implies 
about galaxy outflows in order for predicted gas fractions to be con- 
sistent with the data; further details are presented in Appendix IdI 
How these conclusions are affected by uncertainties in the yield is 
addressed in §|4J2] We then present in §|5]what constraints wind 
metallicity and entrainment fraction considerations place on viable 
outflow models, with a summary and further discussion in §|6] Ap- 
pendix |B] describes the connection between gas masses, accretion, 
and star formation rates in our approach, with implications for star 
formation efficiency. 

Throughout we adopt a cosmology of (H m, ^b, crs, fe) = 
(0.26, 0.047, 0.77, 0.72) and a lChabrieilHoOSah initial mass func- 
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Figur e 1. Mass-metallicity relations listed in Table [T] jKewlev & EllisonI 
120081 and equation [T). The scatter about any given one of these curves is 
0.1-0.15dex, which is much less than the differences in normalization; 
that is, the normalization differences are systematic. The mass-metallicity 
relations in black (T04, solid; D02, dashed) are modeled in detail fi§|4]| 
and Appendix IdI 



tion (IMF), unless otherwise noted. Varying the cosmological 
parameters within t he ranges allowed from observations (e.g., 
iHinshaw et al.l 120091) does not alter our conclusions. The impact 
of varying ilm or 0,b, has, for example, little effect on the shape of 
the M*-Mhaio relation or on the determination of the stellar masses 
in SDSS. Though varying trg does change the number density of 
massive halos, it has little impact on the range of halo masses of 
interest here. Finally, we note that the virial relations only have a 
mild change in normalization when varying cosmological parame- 
ters, without having much impact on our overall results. 



2 RELEVANT OBSERVATIONS 

2.1 The observed z ~ mass-metallicity relation 

Since oxygen is effectively produced only in Type II SNe — the 
deaths of massive, short-lived stars — and H II regions are associ- 
ated with ongoing star formation, the gas-phase "mass-metallicity 
relation" typically refers to only the galaxy's oxygen abundance 
in gas that is currently forming stars; we therefore will use "met- 
als" and "oxygen" interchangeably unless otherwise noted. How- 
ever, though 12 + log(0/H) is measured at the sites of star for- 
mation, the measured abundances are the birth abundances of the 
H II regions; supemovae (the sites of oxygen production) destroy 
their nascent clouds, rendering so-called "self enrichment" of H II 
regions extremely rare. We therefore assume that the galaxy gas 
is well-mixed, i.e., that the mixing time is short relative to the 
timescale for star formation. 

Observationally, oxygen abundance increases with galaxy 
stellar mass. This relation has very little scatter ('^ 0.1 dex in 
12 -I- los[0/H] at fixed stellar mass), though severe outliers do 
exist JPeeples et alJuOOa . l2009f) . The amplitude and slope of the 
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fits to the mass-metallicity relation. 



Table 1. iKewlev & Elhsonl fcOOSi 

where logZg = a + 61ogM* + c(logM*)^ + d(log M*)^, sorted by 
decreasing max(Zg). The two fits we consider in the main text (T04 and 
D02) are in bold. See text for abbreviations. 



mass-metallicity relation, however, are not well constrained, de- 
spite exquisit e and extensive data from the Slo an Digital Sky Sur- 
vey (SDSS; lAdelman-McCarthv et allbOOel) . This ambiguity is 
due to the theoretical uncertanties in how to convert emission-line 
fluxes to 12 + log(0/H), as assumptions must be made about both 
the gas temperature and ionization structure. While the electron 
temperature can be estimated directly using the [O III] A4363 au- 
roral line, this line is extremely weak and usually only detectable 
in very metal-poor environments. Thus, it is common to calibrate 
measurement methods using much stronger forbidden emission 
lines such as [O II] AA3726, 3729, H/3, [O III] AA4959, 5007, 
Ha, and [N II] A6584 based on the so-called direct [O III] A4363 
Te method. However, since [O III] A4363 preferentially emits in 
high-temperature regions, this calibration can lead to an over- 
estimate of the electron temperature base d on this line and thus a n 
under-estimate of the oxygen abundance JKewlev & Ellisoii2008h . 
It is therefore common to instead calibrate strong-line measure- 
ment methods based on theoretical photoionization models. On 
the other hand, there are arguments that such strong-line methods 
over-estimate the true abundance JKennicutt et al.l 120031) . More- 
over, most indicators are either double-valued at low metallicities 
(such as the popular i?23 indicator) or saturate at high metallicites 
as e mission-line cooling shifts to the near-infrared (Bresolin 20061). 
iKewlev & EllisonI ( 120080 highlight many of these issues, 
and derive 12 -I- log(0/H) for a large set of galaxies from 
SDSS using ten indicators (eight of which we consider here: 



Tremonti et al." "2004|^ D02, iDenicolo et"ai 
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Zaritskv et al 



2002 



1994 



KK04, 
KD02, 



Kewlev & Dopital |2002|-M91. iMcGaughl Il99ll ; PP04O3N2 and 



PP04N2, using the IPettini & Pagelll2004i ([ O in]/H/3)/([N II]/H q) 
and [N II]/Hq flux ratios, respectively). The IKewlev & EllisonI fits 
to the mass-metal licity relation ar e giv en in Table | 7] whe re we have 
converted from a lKroupal ( 1200 ll) to a lChabrieJ IJoOSd) IMF and 
from 12 -I- log(0/H) to log Zg, where 



logZg = [12 + log(0/H)] - 12 - log 
= log(0/H)-log 



Mo/Mh 



XMh + YMu. 
15.999/1.0079 



(1) 



0.75 X 1.0079 -K 0.25 x 4.0026 

These mass-metallicity relations are plotted in Figure[T] the scatter 
in Zg at fixed Af* for each mass-metallicity relation is smaller by 
a factor of 2-3 than the spread in normalizations, implying that the 
differences are caused by the systematics discussed above. 

We consider in detail the t wo relations in blac k in Fig- 
urelDandinboldin Table [T] (T04, iTremonti et ai]|2004l and D02, 
IDenicolo et al.ll2002h . The D02 indicator is a linear relation be- 
tween the [N II] A6584/HQ ratio and 12 + log(0/H)calibrated 
against Te metallicities. The relatively low normalization of this 
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Name log M-j 
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K 



f 



7 



8.3298 
8.7265 
9.0892 
9.5141 
9.8941 
10.298 
10.664 
11.053 



0.5153 
0.3084 
0.2062 
-0.07142 
-0.3230 
-0.5548 
-0.8389 
-0.8303 



0.07867 
0.06500 
0.06359 
0.06220 
0.04817 
0.06666 
0.06212 
0.06566 



Table 2. Cold gas fractions logFg = log(Mg/M*) in bins of 
A log Af <- = . 4 dex and the uncert ainty in the mean crios f„ for the 
iMcGauglJ J2005h . iLerov et all <2008ll . and iGarcia-Appadoo et alj J2009l) 
data sets. 



method is common for Te-calibrated indicators. The T04 method 
is based on theoretical stellar population synthesis and photoion- 
ization models combined with a Bayesian analysis of many more 
strong emission lines than used in most methods. While we do 
not favor any one 12 + log(0/H) indicator, we take these two 
mass-metallicity relations as representative of the normalizations 
and slopes observations as a whole. 



2.2 Observed gas fractions of 2 ~ galaxies 

Figure |2] shows how the gas-to-stellar mass ratio Fg (left panel) 
and gas mass Mg (right panel) vary with galaxy stellar mass. The 
open diam onds are total H I gas masses measured from 21 cm 
line fluxes JMcGaughluOOSl) . The crosses a re also H I gas masses, 
with stellar masses measured from SDSS JGarcia-Appadoo et al.l 
l2009l : IWest eTaDl2009l . l2oToh . The filled circles represent the total 
H I + H2 gas masses (including a correction for helium) from The 
H I Nearby Galaxy Survey (THINGS), with the H2 masses derived 
from HERA CO-Line Extragalactic Survey (HERACLES) and the 
Berkeley-Illinois-Maryland Associatio n Survey of Nearb y Galax- 
ies (BIMA SONG) CO measurements l lLerov et alfcOOSh . Though 
there is large scatter in the gas fraction at a fixed stellar mass, gas 
fractions clearly decrease as Ah increases; this be havior is found 
in co somological hydrodynamic simulations (e.g., iHopkins et al.l 
I2OO8I) . The mean logFg in bins of AlogM* = 0.4 dex for 
8.1 ^ logM* $; 11.3 is overplotted with the large solid orange 
squares; we list these means and uncertainties in Table |2] Each of 
these data sets focus on star-forming galaxies similar to those in 
which 12 + log(0/H) is measurable; surveys not restricted to ac- 
tively starfonninggalaxies lead to much lower average gas frac- 
tions JCatinella et al]|201(]|) . 

We parameterize _Fg as power-law of the form 



F^ = 



AL 



AU 



= KfAL, 



(2) 



AL VA'/*,o 

with 7 > 0. Table[3]hsts log Af^.o, Kf and 7 for our adopted gas 
relations. As we show in i)3.1l Fg^ is a more convenient parameteri- 
zation than the commonly used and more arguably intuitive /Xg, the 
gas mass as a fraction of the total baryonic galaxy mass A/* + Afg, 

„ - Mg _ Fg 

^^ - A//g + A/. " 1 + Fg ■ ^^' 

The "total" gas fraction rel ation is a power-law fit t o the combined 
lMcGaughl . lLerov et all and lGarcia-Appadoo et al.l data sets, offset 
by +0.2 dex so that the total gas fractions are greater than those 
implied by the K-S law (see below). In order to understand the con- 



Total 9.6 316228 0.57 

SDSS 6.0 15.85 0.20 

Fiber 2.7 2.24 0.13 

Flat — 0.50 0.00 

Table 3. Gas fraction relation parameters, Fg = A/g/A/* = KfM- 



tribution of a sloped gas fraction relation to the mass-metallicity re- 
lation, we also consider a flat gas relation of Alg — 0.5Af*, shown 
in green in Figure|2] 

For reference, Figure |2] shows how the total baryonic halo 
mass, (f2i,/f2m)Afh, varies with stellar mass (halo mass as a func- 
tion of Aft is calculated as discussed in i]3.2t . The offset between 
the baryonic halo mass and Af* + A'/g is evidence of the so-called 
"missing baryon" problem; the missing baryo ns are either hot o r 
have been expelled from the halos by z = JCrain et al.ll2007h . 
Figure |2] further highlights the fact that for AI^, < 10^° Mq, the 
fraction of baryons in the form of cold gas is roughly constant 
(i.e., the blue and red lines are roughly parallel). Moreover, while 
massive galaxies are gas poor, galaxies with stellar masses below 
~ 10^'^ Mq have most of their mass in the form of gas: the pro- 
cesses responsible for the missing baryons in z — halos must 
also account for this inefficiency of star formation in low mass ha- 
los. We discuss this issue further in Appendix iBl 

The solid line in the right panel of Figure |2] shows the me- 
dian gas fractions obtained by inverting the Kennicutt- Schmidt (K- 
S, Kennicutt 1998; Schmidt 19591) relation, as explained in detail 
in Appendix IaI The shaded contours denote the l-a and 2-a gas 
masses derived for the entire galaxy (Rg — l.l-Rgo,z) in running 
bins of log Af* from log Af* — 8.3 to 11.1; for clarity, galaxies 
falling outside this region are not shown. The solid line is an eyeball 
power-law fit to the median Rg = l.l-Rgo.z gas masses while the 
dashed line is the same for the gas (and stellar) masses within the 
SDSS fiber. The fact that these relations are quite similar to one an- 
other indicates that aperture corrections are relatively small and/or 
that gas fractions are relatively scale-invariant within l.li^go,^. 

The gas masses estimated from the K-S law and the measure- 
ments of total cold gas masses roughly agree with one another 
on the low gas fraction of Fg ~ 0.1 at log Af^. ~ 11, and that 
Fg increases with decreasing stellar mass. The amount of this in- 
crease in gas fraction, however, is in stark disagreement, with a 
range of over an order of magnitude in Fg. The K-S law only 
traces star-forming gas and therefore traces molecular gas more 
clos ely than atomic, a nd dwarf galaxies are deficient in molecular 
gas JLerov et al.l200a) . At large radii in more massive galaxies, the 
gas is predominately atomic, i.e., the H I radii of galaxies is often 
much larger than the op tical (star- forming) radii JBoomsma et al.l 
l2008l ; IWalter et alj|2008h . For the purposes of the mass-metalUcity 
relation, what matters is the total amount of gas that is able to ef- 
fectively mix and dilute metals. A lower limit to this gas mass is 
the gas that is able to collapse and form stars — the gas traced by 
the K-S law. If on the other hand the atomic and molecular gas are 
well mixed (as opposed to, e.g., molecular gas only populating the 
galaxy center and atomic gas being at large radii), then the total 
gas fractions are more applicable. Finally, neither of these gas frac- 
tion estimates include ionized gas; if such gas is not only prevalent 
in typical galaxies but also has efficient mass transfer with both 
supernova ejecta and gas that will cool to form molecular clouds 
(and subsequently H II regions), then even the "total" gas fraction 
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Figure 2. L eft: G as fractions Fg as a function of M*. Rigti t: Gas ma sses Mg as a function of M*. Tlie open black diamonds ai'e H I gas f ractions and 
masse s from lMcG augh (20055); 'he crosses are the same from lWest et al . (2009, 20 ly). The filled circles are H I + H2 gas fractions and masses ILerov et alj 
l200g) . who find that there is very little H2 below log Af, ~ 9.5, which is consistent with the comparison to the H I samples. The red dotted line shows the 
maximum baryonic mass (f2(,/r2,7i)Afhalo, while the green "flat" line shows Alg = 0.5Af*. The blue "total" line is a fit to these data with the normahzation 
increased by 0.2 dex; the orange squares ai'e the mean log Fg of these same data in bins of A log A/* = 0.4 dex with the inner and outer errorbai's denoting 
the uncertainty in and dispersion about the mean, respectively. Gas fractions and masses derived from SDSS data and inverting the K-S law, assuming a radius 
of l.li?9o,z (solid line) and the fiber radius (dashed line); in the right panel, the shaded region corresponds to the 1- and 2-a dispersions in moving bins of 
logA4"*. 



relation will be an underestimate of the gas diluting the galaxies' 
metals. 



2.3 Galaxy outflows 

Though observations of galaxy-scale outflows are notoriously diffi- 
cult, galaxy winds observed in a range of star-forming galaxies dis- 
play a complex, multiphase structur e. Since detectability increases 
with the star formation rate density dVeilleux et al.l2005r) . however, 
the most detailed studies of galaxy winds have been of the out- 
flows associated with extreme starbursts, namely, (ultra)luminous 
infrared galaxies ([UJLI RGs). Studies of blue-shifted absorption- 
lines reveal b oth neutral dHeckman et al.ll2000l;lRuDke et alJuOOa ; 
iMartinI |2005|) and photoionized gas dOrimes et al.l 1 20091) . often 
with several kinematically distinct components. In contrast, X- 
ray emission ai^ound local starbursts such as M82 indicates a hot 
(T ~ 10^-^-10^ K), teniious (n ~ lO-'^-lO"^ ^™~^A '^'"j 
fluid dStricklandfe Stevens! I2OOOI ; IStrickland & Heckmard [20071 
1 20091) . Wind velocities derived from both emissio n and absorp- 
tion line studies ar e typically hundreds of kms^^ dMartinluOOa ; 
iGrimes et alj|2009l) . The outflow velocity v„ of the colder neu- 
tral gas is typically co mparable to on e to a few times the galaxy's 
circular velocity Udrc dMartinl '2005). which is comparable to the 
galaxy's virial velocity iivir (e.g., Diemand et al. 2003). 

The scaling v^ ~ Uvir follows naturally if momentum trans - 
fer from radiation pressure is driving the wind dMurrav et alj2005n . 
For radiation pressure to be effective, the starburst must be Edding- 
ton limited and the outflowing gas has an asymptotic velocity of 



«w(oo) = 2^;csc ( -z 1 ) 

V ^odd / 



1/2 



(4) 



where the escape velocity Uesc is comparable to the virial velocity. 
The wind velocity is there fore typically taken to be v^ = Si^vir. In 
the single-scattering limit dMurrav et al.ll2005l) . 



M^v^ 



istarburst EnucAfsFRC 



(5) 



where Lstarburst is the starburst luminosity and tnuc = 8 x 10~ 
is the nuclear burning efficiency. Thus the mass-loading factoiU -qv. 
is proportional to the inverse of the virial velocity such that 



»7w 



momentum 



MSF 



SOkms" 



(6) 



This same scaling is achieved if the wind is driven by cosmic rays 
dSocrates et alj2008l) . 

On the other hand, the outflow may be driven by energy 
transfer, perhaps from supemovae thermally heating the ISM 
dChevalier & Clegdll985l ; [Pekel & Silklll986| - ISilk & Reeslll998l ; 
iMurrav et alj2005r) . In this popular scenario, 

-M^vl « ^^SN X (7) 

[# of SNe per solar mass of stars formed] Msfr, 



•^ Definitions in the literature of the "mass-loading factor" vary; we take it 
to mean the total outflow mass rate divided by the total star formation rate 
(including short-lived stars). 
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where Es^ ^ \Q erg is the typical energy per supernova and ^ 
is the efficiency with which supemovae transfer energy to the ISM. 
Letting ^ = 0.1, i.e., a 10% efficiency, and taking the number of 
supernovae per unit mass to be 10^^, this yields a mass-loading 
factor of 



energy 



MsP 



73kms" 



(8) 



where we have implicity assumed ii„ ~ Sfvir . While we in general 



for /3 > (or, equivalently. 



consider models in which rj^ 

?7w oc M^^Jj^ , see ii3.2b : it is helpful to keep the normalizations 

suggested by equations ^ and ([Sj in mind. 

Except via the impact of outflows on galaxy gas frac- 
tions (see Appendix |B), the mass-metallicity relation is insen- 
sitive to the total mass outflow rate A'l^. Instead, as we show 
in ij3.1| oxygen depletion due to winds is governed by the 
rate of metal loss, Z„M^, where Z^ is the metallicity of 
the outflow; in our case (see ii2. lb . the mass ratio of oxy- 
gen in the outflowing material. While many metals (oxygen, 
as well as , e.g., iron, s odium, carbon, magnesium, and neon; 
Heckman e t al. 2000; M artini |2005|; IStrick land & Heckman 2007; 



Martin & Bouchai2009i ; .Grimes et al.ll2009, ; .Spoon & Holtii2009D 
are observed in galaxy outflows, there are relatively few observa- 
tions of outflowing oxyg en, and elemental abundances i n the wind 
fluid are rarely reported. IStrickland & HeckmanI ( |2009|) . however, 
find that the X-ray emitting outflow from M82 has a high enough 
metal content that it is consistent with containing nearly all of the 
freshly produced metals in the starburst with an inferred velocity of 
~ 1000-2000 km s^^. Combined with their interpretation that the 
outflow has very little entrained gas (i.e., that it is essentially com- 
prised solely of supernova ejecta), this implies that the metallicity 
of the outflow is quite high. (We note that in this interpretation of 
the data, supernova explosions surprisingly have no radiative en- 
ergy losses when interacting with the ambient ISM [^ = 1 in equa- 
tion [T); see also iHeckmanI |2003| .) This picture is further compli- 
cated by the fact that outflows are likely multi-phase, and the metal- 
licities and escape fractions in, e.g., the cold and ionized phases 
may be different. From the perspective of the mass-metallicity rela- 
tion, however, what matters is the total amount of expelled oxygen 
relative to the total amount of expelled gas, where "expelled" oxy- 
gen or gas is just the oxygen or gas that has either been physically 
ejected from the galaxy or simply heated up such that it cannot ef- 
ficiently transfer mass to the gas that is able to cool and form stars 
and thus be observed contributing to the mass-metallicity relation. 



3 THE FORMALISM 

3.1 The mass-metallicity relation 

The three galaxy masses relevant to the mass-metallicity relation 
are the total galaxy mass in stars, Nh, the galaxy gas mass, A/g, 
and the mass of gas-phase metals, Mz- The model is based on re- 
lating the instantaneous change in these masses via their sources 
and sinks to the instantaneous galaxy star formation rate, AfspR, 
ignoring envir onmental effec t s such as mergers and tidal stripping 
(see also, e.g.,'Tinsley"l974','l980'; Matteucci '20021; iRecchi etall 
2009; Finlator& Dave 200 8; Spitoni et al. 2010). O bservationally, 
Mannucci et al.l <20IOf) and lLara-Lopez et al.l ( 120101) have recently 
shown that Zg has less scatter at fixed Mi, and Msfr than at just 
fixed Af* (i.e., the mass-metallicity relation); there is no evidence 
for evolution of this surface up to jz ~ 2.5. This finding implies that 



the Mi,-MsFK-Zg hypersurface provides a more physical descrip- 
tion of the underlying physics than just the Mi,-Zg^ plane. In the 
formalism, the star formation rate is closely linked with outflow 
efficiencies, and observationally, gas fractions and star formation 
rates are tightly correlated. We review the relevant equations here 
and their derivation in Appendix [Cl 

The mass-metallicity relation is described as 



^w (,a "r 



f'g(i - U 



dlogAfg ^ dlogZg 



dlogAf* dlogA/* 



+ 1 



= 2/ [Cw - Ca + aFg + 1] 



where 



( i ./ rccy ) 



d log A/g , dlogZg 



d log AL d log M^ 



(9) 



(10) 



(11) 



is a factor of order unity. Equation l llOt shows that the metallicity 
Zg is proportional to the nucleosynthetic yield y. Because the IMF 
and Type II supernova yields are high ly uncertain, y is only con- 
strained to be 0.08 < y < 0.023 (e.g.. lFinlator & Davdl2008l) ; we 
adopt a mid-range value ofy = 0.015. We further justify this value 
and discuss models with a varying yield in ii l4.2| 

The denominator of equation JlOt includes terms for metal 
accretion (i^a), metal expulsion (fw), and dilution from gas (aFg). 
The metallicity-weighted mass-loading factors i^a and i^w in equa- 
tion iC6t describe the relative rates at which metals are being ac- 
creted and expelled from the system, and are defined as 



Ca 



ZlGM Afacc 

—^ X 



A'fsFR 

Afw 



ZiGl. 



?7a, and 



Msfr 



(12) 
(13) 



The metallicity of accreting gas, Zigm, is typically taken to 
be zero, though SPH simulations indicate that due to previous 
episodes of enrichment of the intergalactic medium (IGM) from 
metal-cont aining galaxy outflows, the effective Zigm may be non- 
negligible l lFinlator & Davdl2008l ; lOppenheimer et alj|2009l) . Be- 
cause a self-consistent model of an enriched IGM will be based 
on the evolution of the mass-metallicity relation, we will for now 
take Zigm and thus Ca — 0, though we will return to the rami- 
fications of this assumption in §[63] The wind metallicity, Z„, is 
often assumed to be the ISM metallicity l lFinlator & Davg 120081 ; 
lErb 2008) . giving i^w = ?7w. However, Zg is simply a lower-limit 
to the possible outflow metallicity (if the wind is driven by super- 
novae, then it can be metal-enriched relative to the ambient ISM, 
but not metal-depleted). The actual wind metallicity will depend 
on the fraction /o of the outflow that is entrained interstellar gas, 
which has a generic metallicity Zg, and the fraction 1 — /o that 
is from newly e xploded supernovae and th erefore has a metallicity 
^oj,max ~ 0.1 jWooslev & Wcave r|[l995l) . The wind metalhcity is 
thus 



^w — (J- /c)^cj,max ~T~ ./c^g 



(14) 



where we note that /o may vary with galaxy mass and must satisfy 
< /o < 1. 

In this model, as long as galaxies have a given Fg-Ah relation 



(§23), how they got that gas (i.e., ?7a and ?7w) is irrelevant. However, 
for any given wind model r;„, the accretion rate as a function of the 
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star formation rate r^a is uniquely determined. We explore this point 
and its implications in Appendix IbI 

By finding combinations of the yield, outflow strength, and 
gas fractions that combine as stated on the right-hand side of Equa- 
tion dlOt to give Z^{Mi,) on the left-hand side, we can explicitly 
reproduce the observed mass-metallicity relation. This is the tack 
we take in §[4] 



3.2 Connecting galaxies and lialos 

A number of methods have been developed to empirically con- 
nect galaxies to halos. One straightforward approach is the cumu- 
lative matching of galaxy (Wgai) and halo (r t haio) number counts 
JVale & Ostrikej |2004 IShankar et all l2006l : IConrov & Wechslej 
I2OO9I) i.e.. 



ngal(> Mi,) = nhalo(> Mhalo) • 



(15) 



Assuming that each halo (and subhalo) contains a galaxy, equa- 
tion M5\ determines the average mapping between halo mass and 
galaxy mass. 

We a dopt one of the late st determinations of the M*-Mhaio 
relation bv lMoster et alj j2009l top panel of Figure[3ll, 






0.0633(1+2) 



(16) 



Mh 



Mh,o 



-l.06-0.17z 



Mk 



Mh,o 



0.556(l+z)" 



with th e z ero point increase d by 0. 05 dex to correct fro m a lKroupal 
( l200lh to a lChabrierl J2003d) IMF teemardi et alj2010l) . and where 

log Mh,o/M0 = [log 11.88] (1 + z)°-°" . (17) 

The iMoster et alj j2009r) A/^-A/haio mapping is in good agree- 
ment with constraints from galaxy-galaxy lensing, galaxy clus- 
tering, and predicti ons of semi-analyti c models. Following the 
scaling relations in iTonini et alj l l200q , and references therein), 
we have verified that equation ([TTJ yields a TuUy-Fisher relation 
JTullv & F isher 197"?) consistent with the more recent calibrations 
bv lPizagno et al. ( 20071) , as long as the dynamical contribution of 
the dark matter within a few optical radii is less than the one pre- 
dicted by a pure iNavarro. Frenk, & White! ill 996 ) mas s profile, in 
line with many other studies (e.g. .lSalucci et alJuOOTD . Also note 
that the subhalo masses quoted bv lMoster et al.l refer to unstripped 
quantities, which represent more reliable indicators of the intrinsic 
potential well in which satellites formed. 

The iMoster et alj relati on is in br o ad ag reem ent with other 
studies, such as the ones by IGuo et alj (|2010[) and IShankar et alj 
yOOdj), although the latter relied on a stellar mass function based 
on dynamical mass-to-light ratios that cannot be directly used in 
the present study based on SDSS stellar masses. Despite the differ- 
ent techniques adopted, most of the studies find consistent results 
on the A/*-A/haio relation, especially in the mass range of interest 
here, i.e., star-forming galaxies with stellar mass < 2 x lO^^Mn 



and hence halos with mass < 5 x 1 0^^ h ^Mm JFirmani et alj 
l2009l : lMore et alj2010l) . We caution that lNeistein et alj l l201 ll) have 
recently described how the Af*-Mhaio relation could be quite dif- 
ferent from expected by the basic abundance matching technique. 
Allowing satellite galaxy masses to depend both on host subhalo 
mass at infall and on the friends-of-friends (FOF) group mass, 
many distinct galaxy-halo coiTelations are found to satisfy all basic 
statistical and clustering constraints. In particular, successful mod- 




Figure 3. Halo mass, Mh„]„, and virial vel ocity, iJvir. as a function of stel- 
lar mass, Ah,atz = JMoster et al.l2009l) . See §|I21for more details. 



els are found where satellite galaxies are hosted by much lower 
'mass halos than we assume here. 

If winds depend on the potential well depth of the halo rather 
than the mass itself, then the halo virial velocity t^vir is more rele- 
vant than A'fhaio- Roughly speaking, 

GAfhalo 



$, 



(18) 



Rhalo 

where the dependence of the halo radius i?haio on the halo mass 
is a function of both cosmology and the structure of the halo 
JLokas & MamoiJ I2OOII : iFerrareseJ I2OO2I : iLoeb & Peebie3 l2003l: 
iBaes et al.ll2003l) . We connect v^ir to Afhaio via 



fGMi, 



\ ^vir 



112.6 



1/2 



A//h 



i In 



1012 Mq 
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1/3 



(19) 
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1/6 



(1+^: 



1/2 
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where the mean density contrast (relative to the critical density) 
within the virial radius iJvir is A = ISvr^ + 82d - 39d'^, with 
d = n^ - 1, and n^ = n^(l + zf/ \Q ^(1 + z)^ + Oa] 
l lBrvan & Normanll998ljBarkana & Loebl20()lh . The bottom panel 
of Figure [3] shows how v^ir varies with stellar mass in this model. 
We have verified that our Af*-Uvir relation is in good agree- 
ment w ith the Mi,-V2oo relation recently derived bv iDutton et alj 
( l2010al) . 



4 MODELS OF THE MASS-METALLICITY RELATION 

We now turn to what is required to reproduce the observed mass- 
metallicity relation. Rearranging equation l llOt . we find 



Cw - Ca + O-Fg 



(20) 
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where we hereafter take C,a, = (i.e., metal accretion is negligi- 
ble; see ^ 16.31 for a discussion of this choice). Expressed this way, 
the metallicity Z^ is related explicitly to a sum of i^w (a term de- 
scribing outflows) and Fg — A/g/Af* (a term describing the galaxy 
gas content). Equation ( 1201 ) is the principal theoretical result of this 
paper, connecting gas-phase metallicities to gas fractions, outflows, 
and accretion. Functionally, one can use equation l l20b to find work- 
ing models for a given Zg(M^) in several ways: 

(i) Assume y and Zg{Mi,) are known; use trial and error to 
find combinations of (^w(i'vir) and [aFg](A/i) that satisfy equa- 
tion l [20t . 

(ii) Assume y, Zg{A'h), and ("w('yvir) are known; solve for 
d log A^g /d log A'f * in equation ^ and integrate to find Mg ( Af * ) . 

(iii) Assume y, Zg{M),), and Mg{A'h) are known; equa- 
tion J20t then says i^w = y/Zg — 1 — aFg. 

Method [(1)1 works well for developing intuition regarding tensions 
in the data and theoretical wind models, while methods [(ii)| and [(Iir)l 
yield models that exactly produce the observed mass-metallicity re- 
lation. In ii l4. II we explore models with a constant yield y — 0.015, 
focusing in ij |4. 1.2| on what constraining the model to match ob- 
served gas fractions implies about the efficiency of metal expulsion. 
In Appendix IdI we go into some of the more subtle details of how 
different scalings of ^w with Uvir are and are not consistent with 
observed gas fractions. In particular, we show that the no-winds 
model (^w = 0) requires gas fractions that are ~ 0.3 dex higher 
than observed for all galaxy masses, implying that if the yield is 
constant, then the mass-metallicity relation is direct evidence of 
outflows. Finally, in ii l4.2l we consider how variable yields affect 
our conclusions. 



4.1 Models with constant yield 

4.1.1 Choice of Fg- Mi, relation 

Figure|4]shows outflow models (^w(i'vir) for different given Fg-Mi, 
relations [method [(iii)| . As discussed in m.2\ we consider the to- 
tal gas fractions (blue, solid lines), Mg = 0.05(r26/r2m)Afhaio 
(red, dotted), gas fractions as inferred by inverting the Schmidt 
law for SDSS galaxies (purple), and Mg — 0.5A/* (green). The 
A'lg oc A/halo model is included because it might provide a nat- 
ural explanation for the observed turnover in the mass-metallicity 
relation near M* . We find that for the observed normalization of 
Fg(Af^), the slope of the gas fraction relation is largely irrelevant 
in setting the mass-metallicity relation morphology. That is, z — 
galaxies have little enough gas that the mass-metallicity relation is 
shaped by how ("w rather than Fg scales with galaxy mass. This can 
be seen visually in the right-hand panel of Figure |4] at low masses, 
even the flat gas fraction relation has approximately the same (^w 
slope as those models with steep Fg-A/* relations. 



4.1.2 Best-fit models 

We quantify what (^w(wvii) scalings are required in order to re- 
produce the mass-metallicity relation while remaining consistent 
with the observed gas fractions by using method |(ii)[ by taking 
a given (^w we can compare the corresponding Fg to binned gas 
fractions (§|2J2] Table |2]( to calculate a x^ . Parameterizing ("„ as 
(iio/i'vir)"'' + Cw,o, the best-fit model for the T04 mass-metallicity 
relation is ^w = (78kms~^/iivir)~'''*^ + 0.19, as shown in Fig- 
ure[5] We show the A v^ contours for 1-, 2-, and 3-cr using the Ax^- 
to-(7 conversion from lPress et al.l ( Il992n for 5 degrees-of-freedom 
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Cw,0 



T04 


78.0 


3.81 


0.19 


Z94 


63.5 


3.20 


0.23 


KK04 


55.5 


3.04 


0.32 


KD02 


71.0 


3.18 


0.39 


M91 


73.0 


2.47 


0.77 


D02 


79.0 


3.42 


1.25 


PP04O3N2 


90.0 


3.15 


1.50 


PP04N2 


111.8 


2.31 


1.35 



Table 4. Best-fit parameters for ijw = (^^o/^vir)'' + Cw.o the fits to the 
mass-metallicity relation calculated bv lKewlev & Ellison! 120081) and listed 
in Table [T] and the binned gas fractions plotted in Figure [5] These fw are 
plotted next to the con'esponding Zg(Mi,) in Figure|6] 



(8 data points and 3 parameters). The best-fit values do not change 
significantly if the dispersion about the mean is used instead of 
the uncertainties when calculating x'^, and we safely consider the 
points and errors for the binned data to be uncorrelated because 
the measurements for individual galaxies do not depend on one an- 
other We bin Fg instead of Zg because the mass-metallicity rela- 
tion has been more rigorously measured than the Fg-A/* relation. 
The white regions in Figure |5] correspond to models that are un- 
physical because they require negative gas fractions. The best-fit 
models are always close to the border between physical and un- 
physical regions in parameter space, reflecting the fact that gas frac- 
tions at z = are relatively small; it takes only a small change in 
Cw to go from a small Fg to a negative one. 

The best-fit i^„ can be strongly driven by the turnover of 
the mass-metallicity relation and change in slope of the A/*- 
i)vir relation above log Af^. = 10.5. For example, for the T04 
mass-metallicity relation, if we instead only consider the data at 
logA^",^ < 10.5, thebest-fitCw is instead (72kms"V"vir)"*''^ + 
0.41; that is, the velocity normalization vo does not change much, 
but the slope steepens and the constant offset (^„,o increases. 
Whether the best-fit (w shifts to higher b and (w.o (T04 and D02), 
lower b and Cw,o (M91, Z94, PP04O3N2, and PP04N2), or doesn't 
change (KD02 and KK04) when only modeling log A/* < 10.5 
depends on the subtle details of the particular fit to the mass- 
metallicity relation under consideration. In all cases, however, 
Ax^ for the parameters for the best fitting (^„ for a given mass- 
metallicity relation when the entire mass range is modeled fall 
within 1-CT of the log A/* < 10.5 best fitting model for that indica- 
tor (but not necessarily vice-versa, since the best fitting low-mass 
model often requires negative gas fractions if extrapolated above 
IO^'^'^Mq). The l-o- range of Cw for the T04 mass-metallicity rela- 
tion is shown by the shaded yellow and beige regions in the right- 
hand panel of Figure |4] for the log Af* < 10.5 and entire mass 
range, respectively. 

Other metallicity indicators lead to mass-metallicity relations 
that are generally s hallower and have a lower normalization than 
the iTremonti et al.l mass-metallicity relation. This translates into 
i^w + aFg needing to be larger and to scale slightly less steeply 
with mass than seen in Figure [4] the best-fit i^w for all of the 
mass-metallicity relations shown in Figure [T] are plotted in Fig- 
ure[6] (Detailed example models for the shallow, low-normalization 
Ipenicolo et al. (2002) mass-metallicity relation are shown in Fig- 
ure lD2l ) Numerically, observed gas fractions require 2.3 ^ b < 4; 
this scaling with Wvir is much steeper than the canonical models for 
the unweighted mass-loading parameter discussed in ij |2.3l Further- 
more, (^„ must be large (> 1) at all relevant masses. The only way 
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Figure 4. Required fw to reproduce the T04 mass-metallicity relation witli varying gas fraction relations: total (blue, solid), A/g = O.SM* (green, dashed), 
inverting the K-S law from SDSS data (puiple, solid), and A/g = 0.05(f2(,/r2m)A/iiaio (rsd, dotted); see ii l2.2 | for the motivations behind these relations. 
Left: Gas fr actions as a function of stellar mass. The grey triangles in are the gas fractions plotted in Figurel2 HMcGaugtJl2005l : lLerov et alJl2008l : IWest et alj 
I2OO9II2OIOI) and the orange squares are the binned data (§|T2)- Right- log fw as a function of virial velocity corresponding to the gas fractions in the left panel. 
The orange lines are the best-fitting models based on the binned data (see Figure|5); the beige and yellow shaded regions in the right-hand panel show the l-cr 
range in ^w for the entire mass range and log Af* ^ 10.5, respectively. 
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Figure 5. Ax^ contours for the T04 mass-metallicity relation with fw = {'''o/'''vir)~'' + Cw.o- The black "X" marks the parameters with the lowest x^; the 
yellow, green, cyan, and gr ey regions denote solutions with Ax^ ^ l-o", l-cr < Ax^ ^ 2-a, 2-cr < Ax^ 55 3, and Ax'^ > 3-cr, respectively, using the 
Ax^-to-o" conversion from lPress et alj 119921) . The white regions correspond to unphysical (A/g Sj 0) models. 



around a large ^w is if a significant fraction of the gas that is dilut- 
ing the metals is ionized (and thus not included in the observations 
of cold gas our Fg-Mj, relations). 

In the limit of small Fg and large Cw, one c an see from equa- 
tion GOj that Zg oc Cw^ JFinlator & Davell2008l). We are using cu- 
bic fi ts to the mass-metallicity relation (TablefTl lKewlev & EllisonI 
UOOSi) . but for the relevant mass range, the mass-metallicity rela- 
tion has 0.2 < slope < 0.45 for most of the relations plotted in 



relation (Figure [3]l has M* oc v^^^ 
X uiif for log A//* > 10.6). Thus, 



Figure [6] Our M^^-Mhaio-Wv 
for log A'h < 10 (and Af* 
the metallicity Zg is roughly proportional to vl[^ to ii^jj, imply- 
ing that for ("„ ex w^;|?, b should be in the range 1.2 to 2.7. The 
large constant offset (^w,o, however, means that the parameteriza- 
tion presented here (see, e.g.. Figure |5]l cannot be directly inter- 
preted in terms of the simple power-law scalings presented in §[23] 
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Figure 6. Left: The mass-metallicity relation as derived from different metallicity indicators (ii l2.1llKewlev & EIIisonll200a) . relative to the nucleosynthetic 
yield y = 0.015. Right: The corresponding best fitting fw = {'uo/''vir)'' + Cw.o under the requirement that the models' gas fractions are consistent with 
observations. The C,^ parameters are listed in Table|4] 



We also caution that C,vi "^ r]^, and we explore the consequences of 
metallicity-weighting the mass-loading parameter below (§|5](. 

A crucial step in this analysis is the assignment of virial veloc- 
ities to stellar masses. For example. Ipinlator & Dava 1 120081) found 
that C^w oc v~-^ was sufficient to reproduce the z ~ 2.2 mass- 
metallicity relation (which does not differ significantly in slope 
from the shallow relations at z — 0). In their simulations, however, 
Mi, oc Mhaio, which is a shallower relation than our M* oc M^^^^, 
a slope which iMoster et al.l ( 120090 finds to approximately hold to 
z ^ 2 (see their Figure 14). Because Mhaio oc w^ir, these differ- 
ences have extreme consequences for the interpretation of how <^w 
scales with Vvh-- 



4.2 Implications of uncertain or varying yields 

There is increasing evidence that the I MF may va ry with st ar forma- 
tion r at e, and thus with gal a xy ma ss JMeurer et al. 2009; Lee et al.l 
|2009|) . iKroupa & Weidneil | |2003|) suggest that if all stars form 
in clusters and if more massive stars are more likely to form in 
more massive star clusters, the integrated galactic initial mass func- 
tion (IGIMF) depends on the embedded cluster mass function (the 
ECMF, ^cci(Mcci) o c .A^cci , where Afcci is the mass of a cluster). 
iKoppen et al.l ( 120071) showed that for certain choices of ^eci {Med ) , 
the mmax-Afcci relation (where rrimax is the most massive star a 
cluster of mass A/ed can produce), and SN II yields, the mass- 
metallicity relation could be explained without t he need to inv oke 
outflows. Like many previous models, however, iKoppen et alj de- 
rive stellar and gas masses from their star formation rates under 
various assumptions of closed box with inflows. We re-examine 
here the effects of a varying IGIMF on the mass-metallicity relation 
using observed gas fractions. We connect star formation rates and 



stellar masses from observations; the median SSFR of SDSS DR4 
star-forming galaxies can be fit with a power law 



(21) 



log[MsFR/Af*] = -9.83 - O.12(log[Af,/M0] - 10), 

as shown as a histograin in Figure 10 of'Peeples et al. ( 2009| ). 

Like [koppen et al. ( 2007) , we follow Weidner et al.; J2004h 
and take the minimum mass a star cluster can have to be 5M0 and 
the maximum mass to be governed by the galaxy-wide star forma- 
tion rate such that 

log(A4ci, max/Mo) = 4.93 + 0.751og(MsFR/[Mo yr-']). (22) 

We adopt a power-law slope of the ECMF to be /3 = 2 
dLada & LadiJl2003l ; lKoppen et al.ll2007l) . The IGIMF is thus 



Mod, 



ClGIMp('Tl) 



C(m s; mmax(A//eci))Coci(Moci)dAfcci, (23) 



S.OMq 



where ^{m) is the IMF. The oxygen yielqj as a function of 
stellar mass, y{Mi,), will therefore be determined by ^(m), the 
?Timax(Afcci) relation, and the Type II supernova yields, a s shown 
in the top panel of Figure [7] (where we have adopted the iKroupal 
I2OOII IMF). The lKoppen et al.l ( |2007|) yields are sho wn as the solid 
orange line. The purple lines show models with the lLarsonI ( 120030 



<(Afocl), 

. = 1.2Af° 



(24) 



^ The definition of yield used in this section, y = Mox.y /M* , is slightly 
different from the one given in Equation |C7| but that for most purposes 
these are interchangeable. 
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Figure 7. Top: ox ygen yields as a func tion of log A/* for different choices 
of mmax(Meci) iWeldner et alj|2004|. black: p ^arson 200 j, magenta) and 
nucle osynthetic yields jWoosley & Weaveilll995. solid; .Thielemann et alj 
Il99d long-dashed). Effects of different yields on models of the mass- 
metallicity relation as y/Zg = C^w + Cf^g + 1 for the T04 and D02 
mass-metallicity relations are shown in the middle and bottom panels, re- 
spectively. Constant yields are shown in cyan (y = 0.015, solid; ±0.005, 
dotted); gas fraction line types and colors are the same as in previous fig- 
ures, with the total gas fractions as solid blue Hues and Alg = 0.5Mt as 
short-dashed green lines. Details are given in i) l4.2l 



while the black lines show models with the same for mod- 
els with rrama x (Afeci ) derived from the semi-analytic model of 
Weidner etal. (2004, c.f., th e thick solid line in Figure 1 of 
Weidner & Kroupa 20061) ; the IWeidner et al.l relation gives a shal- 
lower dependence of the yield on M^. The thin solid lines show 
models derived using IWooslev & Weaveii ( 11995!) nucleosynthetic 
models (Z = Zpj), while the dotted lines show the same us- 



ing the [Thie lemann et al. ( 1996) models ; as discussed in detail by 
iThomas et al. ( 1998). Thielemann et al.l gives oxygen abundances 
that are higher than lWooslev & Weaveif s. Our fiducial value of y = 
0.015 is shown in cyan; this is ver y similar to the yields from an 
IGIMF with the [Thielemann et al.l models with the I Weidner et al.l 
mmax-Mcci relation. 

The bottom two panels of Figure |7] show how these uncertain- 
ties in the yield translate to uncertainties in outflows when model- 
ing the mass-metallicity relat ion, where we have pl otted y/Zg for 
the various yields and for the lTremonti et al.l ( 120041 . middle panel) 
and lPenicolo et alj 1 120021 bottom) mass-metallicity relations. The 
thick grey lines at y/Zg = 1 denote the boundary below which 
the yields are not high enough to produce the observed metallici- 
ties. As shown in §[3] the gas fractions and outflows must balance 
to give y/Zg, we also show aFg + 1 for two gas fraction models: 
our fiducial "total" gas fractions in blue and a toy Mg = 0.5 A/* 
model in green. The difference between the y/Zg curves and the 
gas fraction curves shows how much outflows are needed. For ex- 
ample, our fiducial yield of 0.015 gives y/Zg that is greater than 
the a/Fg + 1 curves at all galaxy masses, with an decreasing dif- 
ference at high A/*; these differences are what's explicitly plotted 
in Figuresl4land lD2l The dotted y = 0.015 ± 0.05 lines in Figure|7] 
show how these results qualitatively do not change for a large range 
of constant yields; this range roughly shows the uncertainty in the 
yield from uncertainties in the IMF. Note that by exploring a range 
of normalizations of the mass-metallicity relation in tj |4. 11 we are 
exploring a range of y/Zg — the parameter to which our results are 
sensitive. 

The closeness of the orange line (the lKQppen et alj yields) and 
the blue line (total gas fractions) in the middle panel shows that, 
within reasonable uncertainties, outflows are not strictly needed in 
that model. If the normalization of t he mass-meta llicity relation is 
lower, however, then even with the Koppen et al.l yields, outflows 
are required. This can be explicitly seen by comparing the blue and 
orange lines in the bottom panel for the D02 mass-metallicity re- 
lation; moreover, because the difference between these two curves 
is greater at larger galaxy mass, it implies that in this case outflows 
would have to be more efficent at removing metals from massive 
galaxies. Unfortunately, however, the uncertainties in the oxygen 
production of core-collapse supernovae, mmax(Afcci), the IMF, 
and the normalization of the mass-metallicity relation all conspire 
to make drawing strong conclusions over the nature of outflows in 
the case of a varying yield extremely difficult. 



5 OUTFLOW METALLICITY AND ENTRAINMENT 

Supernova-driven galaxy outflows are comprised of some combina- 
tion of supernova ejecta and ambient interstellar medium entrained 
in the outflow. The fraction /c of entrained gas determines wind 
metallicity Z„. As mentioned in §[3TT] the wind metallicity Z™ is 
usually assumed to be equal to the ISM metallicity Zg when mod- 
eling the mass-metallicity relation, but if the outflowing supernova 
ejecta entrains very little gas (which would dilute the wind metal- 
licity) then Zw could be much higher than Zg . 
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We showed in §|4] that models of the observed z = mass- 
metallicity relation are more consistent with observations of 2 = 
galaxy gas fractions when the metallicity-weighted mass-loading 
factor (^w = {Z^/Zg)ri^ scales steeply with the halo virial ve- 
locity, i.e., C^w ~ («o/«vir)'' + Cw,o with 6 > 3. Theoretical 
models for how supernovae drive galaxy-scale outflows, however, 
generally predict that the unweighted mass-loading factor rj^, = 
Mw/MsFR = (ao/vvir) will scale much more shallowly, with 
/3 = 1 or 2 (§|23](. Reconciling these disparate scalings therefore 
requires that Z„/Zg and hence the wind fluid composition varies 
with galaxy mass. 

For any given ("w that reproduces the mass-metallicity rela- 
tion, additionally assuming the form of rj„ (vvir) uniquely con- 
strains the wind metallicity Z^(Ah). If rj„ is constant with galaxy 
mass, then Z^/Zg must increase sharply in lower mass galaxies 
JSpitoni et al.l 1201(1) . Figure [8] shows Zw for the best-fit C^w = 
(TSkms'Vwvir)^'®^ + 0.19 for the T04 mass-metallicity relation 
(left) and Cw = (79kms"V'"vir)^*^ + 1.25 for the D02 mass- 
metallicity relation (right). The dotted, short-dashed, and long- 
dashed lines are for r;^ oc v^-^, v^-^^ , and v^-^^ , models, respectively. 
If 77w has a similar scaling with mass as (^, then Z„ ~ Zg for all 
masses. However, a less steep dependence of r]^ on Uvir implies 
that outflow metallicities should depend less on galaxy mass than 
Zg. Moreover, determining Z„ from galaxy wind observations has 
different systematics than determining r/w, and Zw clearly depends 
sensitively on the scaling of r/w . Figure[8]shows how measurements 
of Zw(A/t) can therefore be used to place unique constraints on 
»7w. 

Physically, different scalings of ?;„ and ("„ (and thus Zg and 
Zw) indicate that the entrainment fraction /e (equation 1 14t varies 
with galaxy mass, offering a clue to the physics of galaxy out- 
flows. If, for example, /o increases with increasing gas mass (and 
thus galaxy mass), it would indicate that the wind fluid does not 
"punch" through a blanketting column density of gas but instead 
sweeps up this material and expels it from the galaxy. On the other 
hand, /c decreasing with increasing galaxy mass, would indicate 
that the ability of supernova ejecta to collect the surrounding ISM 
into the wind fluid depends on the depth of the galaxy potential 
well. We find the former to be the case: to reconcile a steep („ 
scaling with a shallower t^w scaling, then winds driven from deeper 
potential wells must be more efficient at entraining the ambient 
ISM than those driven from shallow potential wells. We also find 
that in order to have the normalization of rj„ be consistent with 
the normalizations suggested in §[23](i.e., vo ~ 70kms~^) then 
the entrainment fraction must be ~ 1, though the exact value is 
dependent on the value of Zcj.max. This is particularly interesting 
in light of interpretations of X-ray emitting outflows in which the 
wind flu id is almost entirely comprise d of supernova ejecta, i.e., 
/e ~ l lStrickland & Heckmanll2009l) . Because iron is primarily 
not made in Type II supernovae and thus likely affected differently 
by star formation driven outflows than a-elements, stellar [a/Fe] 
variations could also be use d to shed light on t he oxygen expulsion 
efficiency of galaxy winds JRecchi et al.ll2009f) . 



6 SUMMARY AND DISCUSSION 

6.1 The approach: modeling a system of galaxies 

We use a simple model of the z — mass-metallicity relation to 
place constraints on star formation driven galaxy outflows. In this 
formalism (§[3T|and equation[9}, the gas phase (oxygen) metallic- 
ity Zg of star forming galaxies is 



Zg = y[Cw-Ca + aFg + i]" 



(25) 



where y is the nucleosynthetic yield, ^a describes accreting met- 
als, ^w describes the efficiency of metal expulsion, Fg describes 
dilution by gas, and a is a factor of order unity (see equation [TTt. 
In the absence of metal accretion ((^a = 0), equation l l25b shows 
that the metallicity Zg is set by a balance of outflows (("w) and gas 
dilution (aFg), with the normalization set by the nucleosynthetic 
yield y. This equation represents a general result: each piece can 
vary with galaxy mass, halo mass, and redshift. To the extent that 
the star formation history is not bursty, i.e., A/sfr varies slowly 
on timescales of lOMyr then the yield y can be taken as constant 
with time, letting equation ( I25l l describe the instantaneous state of 
a sequence of galaxies. Galaxies at 2: = are assumed to live 
on a hypersurface described by their stellar masses, gas fractions, 
metallicities, outflow and host halo properties. By taking gas frac- 
tions and metallicities from observations, we are therefore able to 
uniquely solve for outflow properties in terms of galaxy masses or 
metallicities (that are therefore easily comparable to observations) 
or in terms of the galaxy potential (and therefore easily compara- 
ble to models of the underlying wind physics). The only fitting of 
models to data in this approach is that of funct ional forms to ob- 
servations of the mass-metallicity relation (e.g.. lKewlev & EllisonI 
12008;) and either models or parameterizations to gas fractions as 
a function of stellar mass (ij |2.2b . Because there is theoretical un- 
certainty in which metallicity indicator(s) to use when calculating 
the mass-metallicity relation from data, we do not favor a particu- 
lar indicator when drawing our conclusions, and specifically state 
which constraints come from which pieces of the mass-metallicity 
relation. 



6.2 Resulting constraints 

We consider implications for both the efficiency of star-formation 
driven galaxy outflows and for the content of the outflowing ma- 
terial. The two relevant outflow efficiencies are the efficiency with 
which a galaxy expels its metals, (^ = (Z^ / Zg){AI^ / AIsfr) , 
which we parameterize as (^w = (^o/fvir)'' + Cw,o- The second 
relevant efficiency is that with which a galaxy expels its gas, the 
unweighted mass-loading parameter r/w = M^/Msfr, which we 
similarly parameterize as t]^ — ((To/i'vir) , where /3 is predicted 
to be ~ 1 or ~ 2 with cro = 70-80 km s~^ (§[231(. The content 
of the wind is observed by its metallicity Z^, which can be ex- 
pressed in terms of the fraction of entrained ISM in the outflow, 
/c, where Z^ = {I - /c)Zcj,max + fcZg (equation [14] in §[3lT](. 
Under the assumption that Zigm = 0, we draw the following con- 
clusions by requiring that viable models reproduce both the z — 
mass-metallicity relation and are consistent with observed cold gas 
fractions. 



6.2.1 The necessity of outflows 

Models with no outflows (M„ — ^ (^w ~ 0) are inconsis- 
tent with observed galaxy gas fractions, if the yield y is constant. 
Specifically, in the absence of winds, the gas masses needed to di- 
lute the produced metals are higher at all galaxy masses than the 
total observed cold gas masses; the magnitude of this offset is as 
great as ~ 0.3 dex in Fg = Mg/M*, depending on the particular 
mass-metallicity relation being modeled. 



Constraining star formation driven galaxy winds 13 



-0.5 

-1 

N -1.5 

-2 

-2.5 



o 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 

:t.=([78 km s-i]/vj3,«. + o.i9 : 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 y 

:^„=([79 km s-i]/v,j3.42+i.25 /: 


T04 Zg 


D02 Zg / - 


- Zwi ^w= [40 km s-']/v^.^ - 

: - - ■ z.; ^.=([?o km s-']/v^j^ /: 


- Z^; r,^= [80 km s->]/v,,, / - 


= ---Z.; r7,=([80 km s->]/vj2^ /I 


- - - Zwi ^„=([80 km s-']/v,j3 / - 


: - - Z„; r7^=([80 km s-']/v,,,)? / : 


// 


r / - 


/ ' - 


/ ^ " ,..••■■ - 


y / 


y^-' 




_<::.: - 


*-■""' Zi^"^^ 


.tf.x« 

..- "'^ <^ 


^'' •^'^^ 


- .- ^ - 'X ^ — • 


y.jy^ 


.'^^^.^-^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 



9 10 

log M./M 



11 



o 



9 10 

log M./M 



11 



o 



Figure 8. Wind metallicities Zw for the best-fit fw T04 (left) and D02 (right) mass-metallicity relations (see Figure|6j. The solid fine corresponds to (^w = r\vi 
and therefore Zw = -Zg and /o = 1; different scalings for rjw = (o"o/'''vir) are shown as the dotted (/3 = 1), short-dashed (/3 = 2), and long-dashed 
(/3 = 3) lines. 



6.2.2 Constraints from the normalization of the mass-metallicity 
relation 

Equation l |25l l makes it clear that the nucleosynthetic yield sets the 
normalization of the mass-metallicity relation. From a modeling 
perspective, it is useful to consider the mass-metallicity relation 
normalization relative to the yield (rather than their absolute val- 
ues) because the true nucleosynthetic yield is unknown to a factor 
of two due to unc ertainties in both the IMF and in Type II super- 
nova physics (e.g. jThomas et al.il998i) . Likewise, the overall nor- 
malization of the mass-metallicity relation (i) l2.U is unknown at 
the ~ 0.3 dex level[j In the constant y framework. The normaliza- 
tion of yjZ^ sets the value of the constant offset (^„,o > (which 
is set by the turnover of the mass-metallicity relation, see below). 
The typical required velocity normalization wo ~ 70-80 km s^^ is 
consistent with expectations. 

Low normalization mass-metallicity relations require <^„ > 1 
for all relevant masses; if the true nucleosynthetic yield is larger 
than our fiducial value (y > 0.015), then the efficiency with which 
galaxies expel metals will have to be even stronger. Thus if nor- 
mal quiescently star forming galaxies are not expelling winds with 
Cw 3> 1, then the data prefer a low nucleosynthetic yield and a high 
normalization of the mass-metallicity relation. Furthermore, be- 
cause the mass-metallicity relation shifts to lower normalizations at 



^ Though neither the nucleosynthetic yield nor the normalization of the 
mass-metallicity relation are well determined, the scatt er in log Zg at fixed 
M* is known to be ±0.1 dex JKewlev & EUisorJIlOOSl) . In light of the for- 
malism presented here, this small scatter implies that either the scatter in 
both aFg and fw are small, or they are highly correlated. 



higher redshifts, galaxies at these epochs must have either stronger 
winds or higher gas fractions than their z = counterparts. 

In i) l4.2l we explored the possibility that the nucleosynthetic 
yield y could vary with star formation rate and thus Af* . While the 
possible y(M*) relations are still highly uncertain and the models 
are much more susceptible to the F^-hh relation, we find that a 
wide range of such models have outflows that are more efficient in 
high mass galaxies than in lower mass galaxies. Observations of 
such a trend would be compelling evidence that y varies strongly 
with galaxy mass. 



6.2.3 Constraints from the morphology of the mass-metallicity 
relation 

The morphology of the mass-metallicity relation has two main fea- 
tures: the slope below --^ M* and the turnover at higher masses. 
In the constant y case, the slope of the mass-metallicity relation 
largely determines how C,vi scales with galaxy mass, though with 
some degeneracies with the normalization and constant offset. For 
small Fg, as is the case at z = 0, (^w should scale roughly as 
Z'^^ ~ M^r"'' to Af^r"'*. The power-law scaling of (^„ with re- 
spect to Wvir is typically & ~ 3. The required scalings with respect 
to M^: are fairly robust, while the scaling with respect to i^vir de- 
pends on our assumed Af*-Afhaio relation; if this relation is sig- 
nificantly different from that derived from the abundance matching 
technique, then the C-w <x v^^^ scaling might be able to be relaxed. 
Regardless, this need for a high and mass-depend ent wind effi- 
ciency is in broad agreement with previous s t udies |Dekel & Wool 
I2OO3I ; button et alj|2010bl ; ISawala et al.ll20IOl; ISpitoni et al.ll2010h . 
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The tumovei|j in the mass-metalhcity relation at log A/* ~ 
10.5 may be an observational artifact of the metallicity indica- 
tors saturating at high Zf, (§|2!TJ; however, if oxygen abundances 
do asymptote to a particular value at high masses, then this be- 
havior can be used to place strong constraints on galaxy out- 
flow properties. Specifically, both the normalization of the mass- 
metallicity relation relative to the yield (max[Zg/j/]) and the ef- 
fects of Vv\v increasing sharply above M* ^ lO^^M© (Figure[3j 
must be then taken into consideration; moreover, the interplay be- 
tween these effects can place stronger constraints on viable mod- 
els than just considerations of the mass-metallicity relation below 
10^°'^ Mq. Morphologically, a turnover in the mass-metallicity re- 
lation means that either aFg or C,-w cannot be approximated as a 
power-law. Because cold gas fractions are observed to roughly fol- 
low a power-law with respect to M*, then ^w needs a constant offset 
Cw,o ~ 0.2-1.5, depending on which indicator is used to calculate 
the mass-metallicity relation and/or the yield. In several cases, if 
Mg oc (f2(,/r2m)A/haio> then C,-w can be described as a power-law; 
physically, this would imply that galaxies above M* have large 
reservoirs of ionized gas that are able to efficiently transfer mass 
with colder, star-forming gas. 

If C,vi and 77w scale differently, then the fraction of entrained 
ISM in the wind fluid will vary with galaxy mass. Observationally 
this will be seen as Z^/Z^ varying with mass. As the morphology 
of the mass-metallicity relation constrains the scaling of C,-w with 
mass, the scaling of Z^ and thus t^w with mass therefore depends 
on the slope of the mass-metallicity relation. For example (see Fig- 
ure[8}, for a fixed r/w, a steep mass-metallicity relation will lead to 
a shallower Z^-Mi, relation than a shallower mass-metallicity re- 
lation will. However, since current uncertainties in the slope of the 
mass-metallicity relation are smaller than uncertainties in how (or 
if) 77w scales with mass, measurements of Z„ across a large range 
in galaxy mass, especially above M* , will be particularly useful for 
constraining how t^w (and C,vi) scale. 

6.3 The role of metal-(re)accretion 

Ai z = Q, the assumption that accreting material has a negligi- 
ble metal content (i.e., that Zigm ~ and therefore ^a ~ 0) 
may not be entirely safe. The IGM is enriched as early as z > 3 
( ISongaila&Cowidll996l ; lEUison et al.ll200(]| ; ISchave et alj|2003h . 
and if this material is re-accreted onto galaxies at later epochs it 
could have a significant effect on the shape and normalization of 
the z — mass-metallicity relation. The re-accretion of winds (i.e., 
gas with ZiGM > 0) is a signi ficant component of accre ted gas in 
cosmological SPH simulations JOppenheimer et al.ll2009l) . Though 
the total accretion rate scales with halo mass (Mace oc Afhaio oc 
w^ij, see Appendix IbIi. the contribution of accre ted metals to the 
mass- metallicity relation may not scale so steeply JFinlator & Davg 
l2008h . Moreover, an extra source of metals (^a will imply that the 
overall efficiency of outflows (i.e., the amplitude of (^„) will need 
to be even higher than the ones presented here. However, the reac- 
cretion of wind material seen in SPH simulations may be sensitive 
to numerical issues in the wind implementation; more detailed in- 
vestigations are needed to verify the importance of wind-recycling. 
The metal budget available for re-accretion depends on both the 
amount of metals expelled at higher redshifts and the recyclying 
timescale. We will address the metal content of winds at 2 > as 



^ The turn-"up" at low masses for the Z94 mass-metallicity relation is un- 
physical and due to the cubic fit to the data. 



implied by the evolution of the mass-metallicity relation in a later 
paper. 
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APPENDIX A: INVERTING THE K-S RELATION 



The o bserved Kennicutt- Schmidt (K-S. lKennicutllll998l : ISchmidtl 
Il959r) relation is commonly used to indirectly estimate gas 
masses in star-form ing galaxies in chemical evolution models (e.g., 
ISpitoni et alJuOlOf) , when direct gas masses are expensi ve (or cur- 
rently impossible) to achieve (such as at hi gh redshifts. lErb et al.l 
l2006bf) . or for large samples of galaxies (e.g.. lTremonti et al.l2004h . 
Furthermore, since 12+log(0/H) is measured only in star-forming 
gas, it is reasonable to consider gas fractions that trace this same 
gas. The purple lines in Pigure[2](see also Pigures|4l and lD2l l are 
the gas masses we derive from applying the K-S law to star-forming 
Data Release 4 SDSS galaxies with z - band magnitude errors of 
< 0.01 mag brinchmann et alj IJOOJ : lAdelman-McCarthv et al.l 



120061) . Specifically, we relate the star formation rate surface den- 
sity EsFR to the gas surface density Eg by 



A/sFR 



= A'gE^ 



(Al) 



= 1.67 X 10" 



1.4 



JMqpc-^ 



Mq yr kpc 



fror nlKennicutd (119981 ) , wher e we have corrected for the fact that 
the Brinchmann et al.l j2004l) st ar formatio n rates are based on a 
i Ki'oupa ( 2001) IMP while the [Kennicutj relation is based on a 
ISalpeterl ( 119551) IMP. SDSS spectra are taken within a 3" aperture; 
therefore, to measure total galaxy properties (e.g., star formation 
rates and stellar masses), the fact that the aperture does not sub- 
tend the entire galaxy must be corrected for. We therefore consider 
EsFR and A/* both for the full galaxy-light radius (which we take 
to be 1.1 times the 90th percentile z-band isophotal radius -Rtjo.z) 
and only within the fiber, i.e., we take 



Ag = TlRg — 7r_Riigi-jt = TT X 



The galaxy gas mass is then simply 



190, z, 



solid lines, 
dashed line. 



(A2) 



Ms 



MsFR X 



AS 



Ka 



I/O 



(A3) 



APPENDIX B: OUTFLOWS, INFLOWS, AND STAR 
FORMATION: GETTING THE GAS MASSES 



As shown in SBTT 



MsFK. + Mrccy " M„ 



A/g = A4cc 

= A/sFR(r;a - 1 + /rccy " 7?w) 



and 

dMg 
dAf* 



^a ^w -L I Jrccy 

J- /rccy 



= ^g(l-7). 



(Bl) 
(B2) 

(B3) 



In §|4]we assumed an Fg-A'h relation existed and that as galaxies 
evolve they remain on such a relation. Here we consider, for a given 
r/w, what implications such a relation has on the gas accretion rate 
and how efficiently galaxies are able to turn this accreted gas into 
stars. The above equations imply that the gas inflow and outflow 
rates must be balanced by 



r?a - r?w = (1 - /rocy)-Fg(l ^ ?) - /rccy + 1- 



(B4) 



Thus, for a given combination of r;„ and Fg, we can uniquely de- 
termine rja = A/acc/A/sFR, i.e., the efficiency with which a galaxy 
turns its accreted gas into stars. Por example, if the star forma- 
tion rate is higher than the accretion rate (log rjn < 0), then the 
galaxy is forming stars more quickly than it is accreting gas, i.e. it 
is very efficient at forming stars. We plot log ?7a for the no wind, 
77w = [70kms""'^]/«vir, and ri„ = ([TOkms^-'J/uvir)^ cases as a 
function of stellar mass in the upper-left panel of Pigure lBTI for the 
total gas fraction relation (see Table[3]and Pigure|2ll. Strikingly, r/a 
always decreases significantly with increasing mass — even in the 
absence of winds (solid blue line). This behavior follows directly 
from the steepness of the gas fraction relation (eauation lB4t . When 
outflows that preferentially remove gas from low-mass galaxies 
(r^w oc II ~J, green dashed line; t;~j^, red dotted line) are taken into 
account, r/a likewise increases and steepens to compensate. There- 
fore, while winds may affect how star-formation efficiency varies 
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Figure Bl. Gas accretion rates, star formation rates, and cold gas accretion fractions as a function of stellar mass with varying outflows and specific star 
formation rates. All panels assume the total cold gas fractions described in ii l2.2l Panel (a) shows how r)^ varies according to equation IB4t for different rjw 
models: no wind (solid blue), a momentum-driven scaling (green dashed), and an energy-driven scaling (red dotted). In all cases, high mass galaxies ac crete 
less gas per unit s tar formation than less massive galaxies. Panel (b) shows the expected range in Mace, tot/MspR between the iNeistein et al.1 j2006h and 
iGenel et alj 120081) Mace models (shaded regions) and with three scalings of Mgpjj^/Af* with stellar mass: constant (blue), oc /Xg (green), and the median 
values from SDSS (red). These Mace, tot/Mgpi^ are qualitatively similar at low masses to the r^a shown in panel (a), but increase rapidly at high masses. 
Panels (c) and (d) show the ratio /cold of these two estimates, with varying r^w and the SDSS SSFRs and with varying the SSFR and no winds, respectively. 



as a function of galaxy mass, they are not necessary to explain the 
trend, implying that additional physics is at play. 

This analysis does not entirely reveal what drives the r/a-A/i 
relation. However, the nature of 77a can be unraveled by appeal- 
ing to MsPR and Mace from independent sources. For example, as 
shown in Figure lB2l the median specific star formation rate (SSFR, 
Msfr/M*) in SDSS DR4 star-forming galaxies decreases with in- 
creasing A/i, though there is large scatter in the SSFR at fixed A/*. 



We consider here three scalings for how the SSFR may vary with 
Mi,. The median SSFRs from SDSS are shown as the solid line in 
Figure lB2l A physically-motivated way to have the SSFR to de- 
crease with mass is to postulate that it is proportional to the total 
gas fraction, /^g. The blue dashed line shows /^g x 4 x 10^^ yr 
for the total gas fractions, while the purple dashed line shows 
^g X 1 X 10"® yr for the SDSS gas fractions (note that the SDSS gas 
fractions were derived largely from these same i\/sFR and thus this 
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Figure B2. Specific star formation rates. Tlie shaded regions 1- and 2-(t 
dispersions in running bins o f log Mj, of the aperture- con'ected specific star 
formation rates from SDSS IBrinchmann et al.ll20Q4l) ; the black solid line 
is a power-law fit to median (equation [21). The purple dashed is the SDSS 
/ig X 1 X 10~^ yr and blue dashed line is the total /Xg X 4 X 10~^° yr; these 
offsets imply a star formation timescale of 1-2.5 Gyr. The shaded regions 
are dotted lines are constant J\/sfr. 



is a somewhat degenerate comparison). Finally, we consider a con- 
stant SSFR, A/sfr/A/* = 2 X 10""' yr (log[AfsFR/A/*] = -9.7 

in Figuie lB2t . 

Using extended Press-Schechter theory. iNeistein et al.l |200g) 
parameterize the baryonic accretion rate onto halos by 



A4, 



7.23 



(t^)'"(5Si) <'-'"'"---•<"'' 



where /;, = VlhlQ.^. iGenel et al.l 1200a) find a similar accretion 
rate of dark matter onto halos in the Millineum Simulation, which 
implies a baryonic accretion rate of 

1.07 



A'/a, 



6.34 



A/h 



IQI^Mq 



0.181 



(l+2)'-'M0yr-\(B6) 



Figure B3. Star formation efficiency A/sFR/J\/g as a function of M,,, tak- 
ing A/g to be the total cold gas masses (thick lines, ii l2.2t and Afg = 0.5A/* 
(thin fines) and three choices of the specific star formation rate: constant 
(solid blue lines), oc /.tg (dashed green lines), and the median values from 
SDSS (dotted red lines). In all cases, a steeply decreasing Fg-M^, relation 
is required for the star formation efficiency to increase with stellar mass. 



A/haio/A/* ~ Af*, which is close to the observed rya-Af* slope at 
high masses. This combined double mass-dependent behaviour of 
77a with stellar mass produces the characteristic "U" shape observed 

in Figure IbTI 

The INeistein et al.l and iGenel et al.l estimates of A//acc, tot are 
for baryonic accretion into the halo. However, only a fraction of this 
infalling gas may be usable for star formation; for example, if this 
onfalling gas is shock-heated as it is accreted, then it will neither be 
detected in H I-I-H2 observations nor contribute towards star forma- 
tion (since we are sensitive to J^a rather than A^acc, tot proper, the 
gas participating in star formation is relevant). Therefore, to better 
characterize the fraction of gas that is accreted "cold" — and there- 
fore able to further cool and form stars — we combine the estimates 
of Mace, MsFK, and rja, defining this cold fraction as 



These accretion rates are for matter being accreted into the halo, not 
the galaxy, and can be safely considered as upper limits to Afacc • 

The range of Afacc,tot/A/sFR allowed between these two 
Afacc models and three SSFRs (constant, solid; oc /ig, dashed; 
SDSS median, dotted) are plotted in the top-right panel of Fig- 
ure IB II At low stellar masses, Af* oc Af^^f^ (equation [TT] 
and Figure [3](, which when combined with the nearly linear 
mass-dependence of the accretion rate with halo mass, provides 
Macc/MsFB. ~ Afhaio/A/* ~ AC°-^, which is the appoxi- 
mate trend found at Af^, < 10^° Mq. Equations l [B5l > and ( |B6l l 
state that the overall "efficiency" of mass accretion Afacc /Afhaio is 
roughly constant with halo mass. Therefore, although the host ha- 
los of lower mass galaxies accrete a proportionally equal baryon 
mass, they are less capable at converting this gas into stars. 
At high masses, however, the opposite is true; galaxies become 
more efficient at converting accreted gas into stars. For Af* > 
IO^^Mq, A-L oc ACfo (equation [nj, implying Afacc /A'/sfr ~ 



/cold = Va 



AfsFR 
Aiacc. tot 



(B7) 



where Afacc, tot and A/sfr are generally defined. For illustrative 
purposes, we let Afacc. tot be defined as in equations l iB5t and l iB6i l. 
Note that to be physical, ^ /cold ^ 1- The lower-left panel 
of Figure IBll shows how /com varies with different ?7w scalings, 
assuming the median SDSS SSFRs, while the lower-right panel 
shows how /cold depends on the SSFR in the absence of winds. 

There are several interesting behaviors in the lower panels of 
Figure IbTI worth noting. First, the morphology of /coid(Af*) is 
fairly robust against variations in the SSFR and ri„: it is roughly 
constant, perhaps with a slight rise, for log Ah < 10.5, i.e., below 
about AI* , and then drops precipitously at higher masses. Physi- 
cally, this is a restatement of galaxies with masses near Ad* being 
more efficient at turning gas accreted by their halos into stars, rel- 
ative to either m ore or less massive galaxies IShankar et al.ll2006l ; 
IGuo et al.ll2010h . 
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Second, /coid(A^*) ~ 1 for low-mass galaxies. At face value, 
this would imply that all the accreting gas is available for star for- 
mation. This closely resembles so-called "cold-mode" accretion 
scenario in which gas falling into lower mass halos along filaments 
do not experience significan t shock-heating, thereby easily accret- 
ing onto the central galaxy ([Dekel & BimboimI 120061 : iKeres et al.l 
I2OO5II2OO9I : button et al.ll2010bl) . At higher masses, on the other 
hand, accreting gas may be shock-heated and subsequently unable 
to cool and contribute to star formation. Despite this neat picture, 
however, we find it intriguing that /cold (A/*) is so close to unity at 
low masses. Figure [2] clearly shows that Af* + Afg in these same 
galaxies falls short of accounting for all of the baryons in the halo 
by at least a factor of two. Thus, a large part of the accreted baryons 
must be removed from the halo via strong winds, even if the star 
formation is reasonably inefficient in these galaxies, possibly in- 
duced by a particularly strong supernova feedback efficiency. 

Finally, Fi gure lB3l builds on this analysis to show the star for- 
mation efficiency, traditionally-defined as A'/sFa/A/g, as a function 
of stellar mass for the total cold gas fractions and the three choices 
of SSFR. In all cases, star formation is more efficient in more mas- 
sive galaxies^_they_are forming more stars per unit gas (though 
see ISchiminovich et al. 2010). Several previous analyses of the 
mass-metallicity relation have suggested that a varying star forma- 
tion efficiency with galaxy mass i s required in order to reproduce 
the m ass-metallicity relation (e.g., lBrooks et al.ll2007l : ICalura et al.l 
uOOgh . Figure lB3l shows that this condition is implicitly passed as 
long as gas fractions are decreasing with galaxy mass and star for- 
mation rates vary reasonably with stellar mass, as is observed for 
z = galaxies. We note, however, that with proper choices of ^„, 
the mass-metallicity relation is theoretic-ally able to be reproduced 
with a constant Fg and therefore constant star formation efficiency. 



APPENDIX C: DERIVING THE MASS-METALLICITY 
RELATION 



The instantaneous change in the stellar mass, 

M* = AfsFR - A>/,ccy 

= AfsFR(l-/rc3cy), 



(CI) 
(C2) 



is given by the creation of stars (A^sfr) and the rate at which stars 
return mass to the ISM when they die, A/rccy- (We include the mass 
of stellar remnants in A/*.) The relative rate of these two effects, 
/recy = Afrocy /AfsFR, depends ou the star formation history and 
therefore varies somewhat with time; its effect on our results, how- 



ever, is small, and we are safe to adopt /rt 



0.2. 



The gas mass is also regulated by the star formation rate and 
/rccy, with gas accretion adding gas and outflows removing gas 
from the system. The instantaneous change in Afg is therefore 

Mg = Afacc - AfsFR + Afrecy - Af„ (C3) 

= AfsFR(77a-l + ,/rocy -r?w), (C4) 

where Afacc is the gas accretion rate and Afw is the mass rate of 
outflowing gas. As introduced in i]2.3l we define the mass-loading 
factor 7?w as A/„/A/sfr; analogously, 7?a = A/acc/AfsFR. The 
sources and sinks of metals are essentially the same as for gas, ex- 
cept that each component can have a different metallicity. Hence, 

Mz = ^IGMA>/acc - ^gAfsFR + ^cjA^.ocy - ^wAf„ (C5) 

= AfsFR(y + ^g(Ca-Cw-l)), (C6) 

where Zigm is the metallicity of accreting gas, Zg is the ISM 



metallicity, Zcj is the metallicity of gas being returned to the ISM 
by dying stars, and Z^x is the metallicity of outflowing gas. The 
yield y is the nucleosynthetic yield, which is defined as the rate at 
which metals are being returned to the ISM relative to the current 
star formation rate, i.e.. 



Af, 



y 



new metals 



Afr. 



Zcj /r< 



(C7) 



After the first generation of Type II supernovae explode (~ 10 yr), 
y is constant for continuous star formation; we adopt a mid-range 
value of y = 0.015. 

Since we are interested in the mass-metallicity relation at 
z = 0, and not its rate of change, it is useful to eliminate the time- 
dependence in equations ( ICIMC6I 1. We assume galaxies live on a 
hypersurface of Afg, A'/*, Zg, halo, accretion and wind properties. 
Dividing out the time-dependence in these equations allows us to 
solve for the shape of this surface, with observations setting the 
amplitude. Combining equations l |C2t and l|C4), 



dA/g 
dAfT 



??a - '7^^ 



l+/r< 



/r. 



= Fg{l-^), 



(C8) 



where we include dAfg/dAf* = Fg(l — 7) based on our parame- 
terization of fg — Afg/Af* (equation[2]l introduced in i]2.2| 
The rate of change of the gas phase metallicity Zg is 



d Af z Mz 



-A4 



1 

Mg 



(C9) 



" dt A/g A/g A/s 
We can now combine equations iC2t . l lC4t . l lC6t . and l |C9l l to find 



dZg 
dAf* 



y 



C,a C,v, 



1 



A/g 
A/sf: 



A/g(l - frccy) 
J/ + ^g[Ca-Cw-l-fg(l-7)] 



(CIO) 



(Cll) 



Afg(l-/recy) 

Equation dClOt can be integrated with respect to M-^ to find 
Zg{Ah). Furthermore, using the iKewlev & Ellisonl 1 120081) fits 
(i) l2.1l Table [T), we can turn the problem around: by assuming we 
know the mass-metallicity relation (and dZg/dA'h), we can infer 
the required relation between, e.g., Fg and ^w. Specifically, by re- 



uation lIClOl. we find 






Cw - Ca + 


(C12) 


^g(i -^--^^l^dlogAf, ' dlogAf., 


+ 1 


-1 



= J/Kw-Ca + aFg + 1] ' 
where 
^-n f w'dlogMg dlogZg 



^ d log Mt d log AL 
is a factor of order unity, as given in §[3TT 



APPENDIX D: GENERAL MODELS OF THE 
MASS-METALLICITY RELATION 



(CI 3) 



(C14) 



Following methods [(ii)| and [(lil)l in §|4] Figures |DT] and |4] show dif- 
ferent combinations o f gas fractions a nd ou tflow efficencies that 
explicity repoduce the iTremonti et al.l ( |2004|) mass-metalli c ity re- 
lation; Figure ID2I shows the same for the iDenicolo et al.l 1 120021) 
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Figure Dl. Required gas fractions to reproduce the T04 mass-metallicity relation with varying power-law slopes of CwC'^vir): Cw = (orange, solid), 
[50kms^^]/iivir (pink, long dashed), ([85kms~^]/i)vir)^ (blue, short dashed), ([SSkms^^J/Dvir)"^ (green, dotted); these normalizations are chosen to 
give gas fractions that are as compatible with the observations as possible. Note that all models fit data in (a) and (b) by construction: panel (a) shows the 
T04 mass-metallicity relation (black, solid) and models (colored lines) while panel (b) shows log[(^w + oiFg\ as a function of stellar mass for the four models 
(colored lines) and log[j//Zg — 1] for the T04 mass-metallicity relation in black. Panel (c) shows the model log Fg as a function of stellar mass (colored 
lines) and the observed gas fractions as grey triangles; these are the same observations plotted in Figure[2l jMcGaughl2005l ; lLerov et al.l2008l ; IWest et alj2009l 
I2OI0I) . The model log C,^ as a function of virial velocity are plotted in panel (d) (the i^w = case is unplotted because of the logarithmic fw axis). 



relation. In the top two panels of Figures IDll and ID2I the obser- 
vations are shown as the solid black curves; the colored lines de- 
note models with different scalings of ("w with iivir- Panel (a) shows 
the mass-metallicity relation (log Z^ as a function of log A/*). The 
models are chosen so that they give C,vi + otF^ to equal the ob- 
served y/Zg — 1 (panel b). Gas fractions and C,vi{vvii) are plot- 
ted in panels (c) and (d), respectively. Because of uncertainties in 
the nucleosynthetic yield, the normalization of the mass-metallicity 



relation, and possible saturation of metallicity indicators at high 
12 + log(0/H) (ii l2. It , we will consider both the mass-metallicity 
relation across the mass range 8.1 ^ log M, ^11.3 and restricted 
to below A/* ~ 10.5 Mq. 

The gas fractions needed to dilute the metals in the absence of 
winds (("w = 0) are shown as the solid orange line; these gas frac- 
tions are higher by a factor of > 3 than observed cold gas fractions 
in typical z ~ galaxies. For a non- varying yield, outflows are 
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Figure D2. Same as Figures lDT] andl4l but for the lDenicolo et al.1 J2002h mass-metallicity relation. The normalizations for (^w ex: t'~ ^ in the middle two panels 
are chosen to give gas fractions that areas compatible with the observations as possible and are: [185kms~^]/i)vir (pink, long dashed), ([100kms~^]/Dvir)^ 
(blue, short dashed), ([90kms^^]/t)vir)'^ (green, dotted). 



22 Peeples & Shankar 

therefore required in order to Iceep thie observed mass-metallicity 
relation consistent with galaxy gas fraction observations. This con- 
clusion holds even more strongly for the other mass-metallicity re- 
lations plotted in Figure[T] in the absence of winds, lower metallic- 
ities imply higher gas fractions. 

The other colored lines show the required gas fractions if 
we assume C,vj oc v~^ (pink, long-dashed), oc v~.^ (blue, short- 
dashed), or oc v~^^ (green, dotted). For the T04 mass-metallicity 
relation, both the momentum-driven and energy-driven C,vi scalings 
require Fg to scale more steeply with mass than is observed; lower 
normalizations of (,,„ force Fg to asymptote to the no winds case. A 
steeper scaling of ("w with Uvir, however, leads to more reasonable 
gas fractions. 



